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We consider QCD with two degenerate flavors of light quarks(up and down) at asymptotically high 
isospin (/r/) with zero baryon chemical potential (/rs) and calculate for the first time a quantitative 
expression for the critical temperature of the deconfinement transition in this regime. At high isospin 
chemical potential and sufficiently low temperatures this theory becomes equivalent to a pure Yang- 
Mills theory and accordingly has a first order deconfinement phase transition. Although this was 
conjectured in a seminal paper by Son and Stephanov in the year 2001, the critical temperature 
of this deconfinement phase transition was not computed. This paper computes the energy scale 
associated with this transition as a function of the chemical potential /ij by relating the parameters 
of the equivalent Yang-Mills theory to those of the underlying theory. We also relate the equation of 
state in one strongly interacting regime of QCD namely at finite isospin density to that in pure Yang- 
Mills, with the latter being amenable to straightforward numerical calculation. Our results for the 
critical temperature of deconfinement transition can be compared with future lattice calculations. 


I. INTRODUCTION 

One of the goals of modern nuclear physics is to un¬ 
derstand the rich structure of the QCD phase diagram at 
finite density and temperature. This is because various 
such phases of matter can be realized in nature under dif¬ 
ferent circumstances. The core of neutron stars and the 
evolution of the early universe are two of the most promi¬ 
nent examples where a knowledge of the phase diagram 
would help our understanding. 

It is very difficult to get a handle on this problem start¬ 
ing directly from QCD. Accordingly it is useful to learn 
what we can from regimes which are tractable even when 
they are not of direct phenomenological relevance. In 
this paper we focus on the low temperature regime of 
two-flavor QCD with zero baryon chemical potential and 
very large isospin chemical potential. The study of this 
regime was originated by Son and Stephanov [1, 2] who 
noted that in the limit of extreme isospin chemical poten¬ 
tials this theory became equivalent to a pure Yang-Mills 
theory. This implies among other things that unlike at 
zero chemical potential the theory has a first-order decon- 
fining transition. Although, the first order deconfinement 
transition at high isospin chemical potential was conjec¬ 
tured by Son and Stephanov, the critical temperature for 
this deconfinement transition has not been calculated be¬ 
fore. This paper, relates the properties of this emergent 
Yang-Mills theory to the parameters of the underlying 
theory and calculates the critical temperature of this first 
order deconfinement transition. 

One exciting feature of this analysis is that we can 
relate the equation of state at large isospin chemical 
potential and low temperature to the equation of state 
of pure Yang-Mills theory. Since the equation of state of 
the pure Yang-Mills theory is computable using numeri¬ 
cal studies on a lattice with relatively modest resources, 
this means that the equation of state of QCD at large 
isospin chemical potential and low temperatures is 
effectively computable with modest resources—including 


a computation of the critical temperature of the first 
order transition. 

To put this problem in context, it is useful to re¬ 
call that considerable progress has already been made in 
mapping the phase diagram both in isospin density vs 
temperature plane and baryon density vs temperature 
plane. The biggest hurdle impeding progress in this 
effort is the strength of the QCD coupling constant 
at moderate energy scales. A significant portion of 
the phase diagram which corresponds to observable 
phenomena such as the dynamics at the core of the 
neutron stars or the heavy ion collisions is inaccessible 
to perturbation theory calculations as the theory is 
strongly coupled in this regime. In order to circumvent 
this difficulty numerical calculations using lattice have 
been employed and this has been remarkably successful 
in some cases. One such success story was in cosmology, 
where it answered questions relating to whether the 
universe went through a phase transition while passing 
from a quark gluon plasma phase to a hadronic phase at 
temperatures around 150 MeV. This problem was satis¬ 
factorily handled by finite temperature lattice QCD and 
it was concluded that the universe did not go through a 
QCD phase transition but a crossover as it cooled [3, 4]. 
But lattice QCD fails at finite baryon chemical potential 
at low temperature. The reason behind this is that 
the fermion determinant in the presence of a chemical 
potential for the quark number becomes complex, 
which causes lattice algorithms involving the method of 
important sampling to break down. This is known as 
the sign problem [5-8]. But an understanding of this 
regime of the phase diagram is absolutely necessary for 
neutron star physics. This is because the mass-radius 
relation and the transport properties of a neutron star 
are determined by the phase of matter at the core 
where baryon density is high and the temperatures 
are low [9-13]. Efforts have been made to explore the 
finite baryon density regime using Nambu-Jona-Lasinio 
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model [14-18], Quark-Meson model [19], Polyakov loop 
extended Nambu-Jona-Lasinio model [20-24], Polyakov 
loop extended quark-meson model [25, 26]. But these 
models are not very rigorous. Since the temperatures of 
interest are very low, and degrees of freedom involved 
are fermionic, Fermi liquid theory of nucleons, which is 
a theory of nucleon quasi particle excitaions about the 
Fermi surface is used to make predictions for baryon 
densities that correspond to a Fermi energy of 100 MeV. 
Very high baryon density and low temperature regime 
on the other hand can be dealt with analytically. In 
this regime, it was conjectured that the quarks and the 
gluons became deconfined and the quarks at the Fermi 
surface were weakly coupled which made rigorous per¬ 
turbative calculations possible [27]. The ground state of 
matter at highest baryon densities and low temperatures 
was found to be in a color-superconducting phase known 
as the color-flavor-locked phase where a condensate of 
Cooper pairs near the Fermi surface generated the color 
Meissner effect shielding the gluons. 

Apart from the regimes discussed in the previous 
paragraph, there exists another regime of the phase 
diagram, given by finite isospin chemical potential but 
zero baryon chemical potential. This paper explores the 
confinement dynamics in this regime. There are two 
reasons why this regime is of interest to us. The first 
reason is that in this regime lattice calculations are not 
plagued by the sign problem and hence may be practical. 
This means that lattice simulations can be used to 
test the results of model calculations and extrapolated 
results from rigorous perturbative calculations. In our 
paper we look at asymptotically high isospin chemical 
potential where we use perturbation theory to make 
predictions. Extrapolation of our results to moderate 
isospin densities, where perturbation theory breaks down 
can be compared with lattice results to determine the 
accuracy and the range of validity of the extrapolations. 
The second reason is that this regime may help give 
some insight into neutron stars as the imposition of 
charge neutrality in such objects(neutron star cores) 
leads to the suppression of proton fraction compared 
to the neutrons. This corresponds to a finite isospin 
chemical potential. Despite this, our system differs from 
the realistic scenario in these objects as there is also a 
finite baryon chemical potential in them as mentioned 
before. 

The intention of this paper is to look at the finite 
isospin regime for zero hb ■ This is largely because 
the intractability of a finite hb region using lattice 
simulations imply that our predictions for such a 
scenario cannot be tested for moderate densities. Also 
the introduction of a finite /is in the problem will 
further undermine the effect of the isospin chemical 
potential: the effect in which we are primarily interested. 

Note that our system is unstable with respect to 
weak decay and is not electrically neutral. Hence the 


thermodynamic limit of such a system does not exist. 
But we can ignore the electromagnetic and weak effects 
in order to isolate the strong dynamics [2], There have 
been considerable recent developments in this region of 
the phase diagram [1, 2, 28-34]. 

The phase structure of QCD at finite isospin chem¬ 
ical potential and zero baryon chemical potential with 
two light quarks was first considered in [1, 2], For small 
enough hi, i.e. hi is much smaller than typical hadronic 
scales, chiral perturbation theory can be used to predict 
that charged pions(7r _ ) start condensing for hi > m n 
where ra, is the mass of the pions. On the other hand, 
the opposite limit, where hi is much larger than typ¬ 
ical hadronic scales and anti-up and down quarks form 
Fermi spheres of radius respectively, interactions can 

be described using perturbative QCD due to asymptotic 
freedom. Just as in BCS in ordinary metals, if there is 
an instability at the Fermi surface caused by attractive 
interaction in some color, flavor and angular momentum 
channel, there will be a condensation of Cooper pairs in 
the problem at hand as well. In this high isospin regime, 
it is believed that, due to the strength of the attractive 
perturbative one gluon exchange interaction in the color 
singlet channel that a color singlet condensate of quark- 
antiquark Cooper pairs indeed does form. This means 
that two possible forms of the condensate are and 
ud. The instanton induced interaction favours the con¬ 
densation in ii'ysd channel over the ud one and hence the 
condensate has pion quantum numbers. It was conjec¬ 
tured that there was a smooth crossover between pion 
superfluid at low isospin to quark-antiquark condensate 
at high isospin density [2]. For an asymptotically high 
Hi the phase diagram in /ij — T plane was argued to have 
a first order deconfinement phase transition apart from 
the second order superconducting phase transition at a 
higher temperature. 

The argument is as follows. The condensate at high 
isospin density at temperatures smaller than the gap, 
has pion quantum numbers and is color neutral. This 
ensures that the gluons are not screened by the conden¬ 
sate. For temperatures well below the gap it is also ener¬ 
getically expensive to excite a quark quasi-particle out of 
the Fermi sphere and hence there is no Debye screening of 
the gluons either and we are left with a pure glue theory. 
A pure glue theory with SU( 3) gauge group is known 
[35] to have a first order deconfinement phase transition 
with increasing temperature. Therefore we can expect a 
similar behavior in our problem. It was argued by Son 
and Stephanov that the critical temperature of this de- 
confinement transition should be lower than the critical 
temperature of the superconducting transition [1]. Inter¬ 
estingly, as there is no phase transition with increasing 
temperature for Hi = 0, it was argued in [2] that the first 
order phase transition line has to end at a critical point 
if one assumes quark-hadron continuity along the T = 0 
isospin axis. 

As mentioned above, It is expected that for A > T 
and hi ^ m P the static color charges are not screened 
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due to the absence of both Debye and Meissner masses 

and the theory describing the low energy dynamics of the 

system is governed by the Lagrangian density ££ = 

where g is the coupling constant that matches that of the 

original theory at the scale A. But this is not a complete 

picture. Although the quarks are bound in SU (3) singlet 

Cooper pairs, they can still partially screen the static 

color charges by changing the chromodielectric (dielectric 

constant for the color electric field) constant of the theory 

e from unity. As we will see, in this theory e > 1 and the 

coulomb potential between static color charges is smaller 
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than ^ by a factor of e where r is the distance between 
the static charges. This amounts to having a pure glue 
theory where the speed of gluon is given by < f 
where, A stands for the chromomagnetic permeability. 

While this general picture was introduced some time 
ago [2] to the best of our knowledge, quantitative pre¬ 
dictions of properties of the effective pure glue theory in 
terms of the underlying theory are lacking. A principal 
purpose of this paper is to compute this chromodielec¬ 
tric constant (and the chromomagnetic permeability) as 
a function of gj integrating out the quarks around the 
Fermi surface in an energy range of the order of the gap. 
For asymptotically high /xj, we find that eA 1 where 
e and A have been defined with respect to eo and Ao of 
the vacuum by absorbing the vacuum polarization effect 
in g. We also predict the effective confinement scale as 
a function of /i/. The confinement scale is related to the 
chromodielectric constant in a simple way which is also 
elaborated in this paper. The scaling of the confinement 
scale with gj obtained in this paper should be compared 
with scaling results obtained from lattice simulations if 
possible in future. 

A particularly interesting fact is that there is a first- 
order deconfinement transition at large /i/. We cannot 
analytically compute its critical temperature from first 
principles. However,if we can determine the ratio of the 
deconfinement transition temperature to the scale of the 
the theory for pure Yang-Mills, we can then predict the 
phase transition temperature for large gi. Moreover, as 
noted artier, the calculation of the deconfinement tran¬ 
sition temperature for pure Yang-Mills is computable on 
the lattice with quite modest resources, so that a pre¬ 
diction of the QCD deconfinement phase transition at 
large gi is viable as is a more general calculation of the 
equation of state. The equation of state turns out to be 
different from the equation of state of pure Yang Mills 
theory suitably rescaled. Apart from the low lying glu- 
oynamics, the theory has a meson Goldstone mode which 
contributes to the equation of state. While the mode is 
massless, and thus it contributes at low energies, the scale 
of its interactions is well above scale of the gluodynam- 
ics. Thus, while the Goldstone mode contributes to the 
pressure, it is also effectively non-interacting and does 
not affect the gluopynamics. 

The story of a pure glue theory anisotropic in space- 
time at sufficiently high isospin chemical potential is rem¬ 


iniscent of the ‘two flavor color-superconductor’ (2 SC) 
phase at low temperature and high gs [36]. In case of 
the 2 SC phase at finite baryon density, two quarks form a 
Cooper pair giving rise to a condensate that breaks color 
517(3) down to color 517(2). Five of the eight gluons 
become massive due to the color-Meissner effect of the 
condensate, but the remaining three remain massless. At 
temperatures lower compared to the gap Debye screen¬ 
ing for these three gluons is absent giving rise to pure 
517(2) gluodynamics. However, just like in our problem, 
the quark-loops alter the chromodielectric constant of the 
system. 

The paper is organized as follows: Sec. 2 discusses the 
microscopic Lagrangian that is used to determine the low 
energy constants of the gluodynamics. Sec. Ill outlines 
the evaluation of the chromodielectric constant and the 
chromomagnetic permeability using the microscopic La¬ 
grangian. Sec. IV estimates the energy scale of the de- 
confinement transition. Sec. V discusses the connection 
of the equation of state of QCD at large isospin chemical 
potential and low temperature to that of pure Yang Mills 
and computes the critical temperature for the deconfine¬ 
ment transition followed by a concluding section. 


II. MICROSCOPIC LAGRANGIAN 

In order to quantify the effect of the paired quarks on 
the chromodielectric constant of the system we need to 
be able to write down the effective theory for the glu¬ 
ons. The low energy constants of this theory need to 
be extracted by integrating out quark loops from a mi¬ 
croscopic Lagrangian, which in this case is the QCD La¬ 
grangian. We emphasize again that since we are working 
at asymptotically high densities it is justified to make 
use of perturbation theory. We start with the QCD La¬ 
grangian with two light (massless) quarks in the funda¬ 
mental representation of the SU (3) color group with an 
isospin chemical potential p,/, 

% = + na 0 T 3 )ip - ^ (i) 

where, 

, F“ v = d^-duA^+gfahcA^Al and = d^-igA^U. 
The generators of color 5/7(3) are denoted by t a and r 3 
is the third Pauli matrix given by 


u and d stand for up and down quarks respectively. As 
mentioned above the condensate has pion quantum num¬ 
bers and is of the form uj^d. In order to analyse this 
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condensate we define 



( 4 ) 


where d = 75 d. We intend to find the propagator for the 
ip quarks in the presence of this condensate and to do so 
we introduce an auxiliary field A such that the inverse 
propagator for the quarks is given by 


G~ l 


(j^ + nn° a \ 
l A 7 "Pm - M/7°y 


( 5 ) 


with A = yoAiyo- The quarks are expected to acquire 
a gap in their spectrum which is to say that it is ener¬ 
getically expensive to excite a quark bound in a Cooper 
pair. This gap is characterized by the magnitude of the 
auxiliary field when the action is minimized with respect 
to the auxiliary field. This amounts to solving for the 
gap, A using Schwinger-Dyson (S-D) equation. We do 
not need to solve for A for the purpose of our paper, 
althouh it is straightforward to do so. The gap was esti¬ 
mated in [2] to be of the form A = b\gj\g~ 5 exp(— c/g) 

o 2 

where c = and g was evaluated at the scale gi. Note 
that, in ordinary superconductivity with a four-fermi 
interaction the gap is of the form A ~ exp(— k/g 2 ) 
where k is some numerical constant. Here, however, 
A ~ exp {—c/g). The source of this behavior is that 
instead of a contact four-fermi interaction we have long 
range magnetic interaction, just as in the case of a 
superconducting gap at large gs- b was estimated to be 
~ 10 4 [2]. It is easy to see from the expression for the 
gap that as g gets smaller the gap A becomes smaller 
compared to gj. This is a standard feature of any BCS 
calculation and is the only consistent solution of the 
gap equation other than A = 0. In our problem as gi 
becomes large, the strength of the coupling becomes 
smaller and we obtain gj A. 


More over, it can be expected that at asymptotically 
high isospin densities, chiral symmetry is restored 
asymptotically. The evidence of this comes from [2], 
There it is shown that for gi > m^, the value of the 
chiral condensate is given by < uu + dd >= 2ipip cos(a) 

and that of the charged pion condensate is given by 
_ 2 
< uj 5 d > +h.c. = 2ipipsin(a), where, cost a) = 

There we notice that as the isospin chemical potential 
reaches from zero and grows further, the condensate 
starts rotating from a chirally broken qq phase to a chi- 
rally symmetric charged pion phase. The condensate is 
expected to crossover to a completely chirally symmetric 
phase at high g,j. The imposition of chiral symmetry 
on the condensate gets rid of and 7 from the 

basis. We can finally express our order parameter in the 
reference frame of the medium i.e. u = ( 1 , 0 ) as 

A = iA 4 + ?A 5 7 .fc 7 o ( 6 ) 

= * (A + Ap + A~Ap) (7) 

. Ap and A~ are free quark on-shell projectors given by 

A p = \ (! ± 7o7-P) ( 8 ) 

, where A + and A~ are real numbers. If we restrict 
ourselves to the 1 S'o channel, we have A + = A~. For 
the purpose of this paper we are only interested in the 
1 S’o channel and we do take A + = A~ = A. With this 
structure of the condensate in mind, we return to deriving 
the propagator for the %p quarks. We do so by going to 
Euclidean time which corresponds to a finite temperature 
calculation and then take the limit T —> 0 at the end. The 
inverse quark-propagator in Euclidean action is given by 

~ 1 _ - 7V + M7° A \ 

^ A -iyV - 7 V - At7 o1 / j 

. Inverting G~ 1 we arrive at the propagator 


For the sake of clarity we explore the Dirac struc¬ 
ture of the condensate here to some extent. Forgetting 
for a moment that the condensate is a constant in 
space-time, we write it as a function of momentum as 
A(fc). The gap matrix A(fc) is diagonal in color space 
but not in Dirac space. Just like any 4x4 matrix in 
the Dirac space it can be expanded in the basis of 16 
linearly independent 4x4 matrices. We however, are 
interested in the channel with a total angular momentum 
of zero. If we restrict ourselves to the J = 0 channel, 
the basis includes only 8 Lorentz invariant matrices, 

1 , 7 ^, 7 /^ M , 7 ^ 7 ^, 75 , 757 /^ M , 757 ^, 757 ^ 7 ^" 
where is the 4 velocity of the medium. To restrict 
ourselves more, we are interested only in the odd parity 
channel. This is because, as mentioned before, although 
perturbative one gluon exchange gives rise to equally 
strong attractive interaction in both ud (even parity) 
and W 75 G? (odd parity) channels, the instanton induced 
interaction favors the uj^d channel. The requirement of 
odd parity shrinks the basis to 1 , 


where 


G = 



£- 

G~ 


n- - ST' ~ ik o + (M “ ek ) 

e=± \ e k) 


£"(*) = *£ 

e=± 

s + (fc) = -i£ 


A e A e 

- k2 o ^ (e e fc ) 2 
A e A -e 
- k2 o - « t ) 2 


( 10 ) 


( 11 ) 


and (e |) 2 = (|k| — eg) 2 + (A 6 ) 2 where e = ±. Note 
that, as expected, the dispersion relation of the quarks 
have become gapped. Also, the dispersion relations with 
e = — 1 correspond to anti-up and down quarks where 
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as the ones with +1 correspond to up and anti-down 
quarks. Apart from the propagator of the quarks in this 
basis, that is ip quarks, we also need the interactions of 
the ip quarks with the gluons to make any progress. The 
interaction term between the quarks and the gluons in the 
QCD Lagrangian, i. e. ip / y IJ, A l pt a ip needs to be rewritten 
in the basis of the ip quarks. For the sake of clarity we 
write down explicitly the interaction term between the 
gauge field and ip quarks and it is given by 


7 (-igi^A^t a 

H 0 " 


0 

—ig r ) fl A a t < 


iP 


( 12 ) 


III. EFFECTIVE ACTION 

Armed with the propagator for the quark and the in¬ 
teraction of the quarks with the gluons we can now write 
down the form of the low energy effective Lagrangian 
from symmetry arguments alone. The microscopic La¬ 
grangian had SU[A)color x SU(2)fi avor x U(l) em internal 
symmetries, none of which are broken by the condensate. 
This means that the low energy effective theory should 
respect these internal symmetries. The microscopic La¬ 
grangian however, was not Lorentz invariant due to the 
presence of the chemical potential which broke Lorentz 
symmetry to rotational invariance explicitly. This means 
that the effective Lagrangian for the gluons will be rota- 
tionally invariant but not Lorentz invariant. Apart from 
this, by requiring locality the form of the effective action 
can be obtained as 

S g iue =J^g2 J ^ ( 2 E » 2A Z, “) (* + ° 

, where E a and B a are chromoelectric and chromomag¬ 
net ic fields respectively, g is the QCD coupling at energy 
scale gi. For gluon momentum q , smaller than the gap 
A, higher dimensional operators which include powers 
of q or more powers of E and B fields are supressed by 
powers of ^ and the effective action behaves like a pure 
Yang-Mills theory, e and A are the low energy constants, 
to be identified as the chromodiclectric constant and the 
chromomagnetic permeability of the medium. These can 
be explicitly calculated by computing the polarization 
tensor for the gluons, which amounts to integrating out 
the quark quasiparticles upto energies of the order of the 
gap around the Fermi surface. In other words, we use 
the one loop polarization tensor to obtain the part of the 
effective action that is quadratic in the gluon fields and 
also contains the effects of the integrated out quarks. In 
this action for the gluons, we rewrite the gluon fields in 
terms of the chromoelectric and chromomagnetic fields 
using the definitions 


rpa _ TAd 

(14) 

'k = Gj k Eij 

(15) 


and compare with Eq. 13. The two Lagrangians are 
compared in order to write down e and A in terms of the 
various components of the polarization. Upto leading 
order in the expansion of gluon momentum over the gap, 
at zero temperature this matching calculation yields the 
following relations: 

n“( 90 ,q) = -(e-l)|q| 2 <U 

n^x = - ( e-1)9 ‘> (16) 

n^o,q) - jqpP) = Q - i) \<l\ 2 $abS jk 

where, II is the polarization of the gluons. The indices 
in superscripts on II stand for space-time indices and the 
subscript for color indices. 

We now have to calculate the polarization tensor. At one 
loop it is given by 

n ab(P) =9 2 t Y 1 j 7WW Trace [l l taG + {k)^t b G + {k-p) 

k 0 J ^ 7I ' 

+7 ’t a G~ {k)^ 3 t b G~ (k - p) 
+Yt a '£~{k)'y J t b '£ + (k -p) 
+Yt a 'S + (k)'y 3 t b T,~(k - p)\ 

(17) 


. The expressions for G + ,G ,_ ,E + and E _ are substi¬ 
tuted from eq. 11 in eq. 17. It is easy to see after the 
substitution that 11 ™ = n™ = II 33 and so on. 

At this point we should note that the trace in eq. 17 
is a standard one and especially similar to the unbro¬ 
ken SU( 2) gluon polarization in the two flavor color¬ 
superconducting phase. The reason behind this simi¬ 
larity has been explained before and is related to the 
fact that the low energy dynamics of both the phases 
are driven by gluodynamics. The above mentioned trace 
was calculated by [37] and the integral over 3 momen¬ 
tum was computed in [36] in the context of the 2 SC 
phase. We should remember that the three momentum 
integral is dominated by anti-up and down quarks with 
momenta close to the Fermi surface. Hence, the up and 
anti-down quark contribution to the above integral is ne¬ 
glected which amounts to neglecting terms with denomi¬ 
nators (|k| — e/r ) 2 + (A e ) 2 or (|(k— p)| — eg) 2 + (A e ) 2 with 
e = — 1. Because the trace and the integral in eq. 17 are 
standard, instead of going into the details of the calcu¬ 
lation we just state the results here. We find that as 


T — » 0, around <? 0 ,q = 0 up to leading order 7^72 

mansion n 00 — 9 2 ^J\^\~ j rr n _ 9 2 sjijx 

pansion, n aa — i 8 lr 2 (A +) 2 ana ii ab — i 8 tt 2 (A+) 2 ° ° ab • 
The appearance of p,/ in the numerator is due to the fact 
that the integral in eq. 17 is dominated by momenta 
near the Fermi surface. This means that the three mo¬ 
mentum integral essentially turns into a one dimensional 
integral multiplied by the area of the Fermi surface which 
is proportional to the square of the Fermi momentum or 
the chemical potential. Note that n* J 6 (— S 3k — ^ 5 -) = 0. 


ex- 
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Comparing eq. 16 with these expressions for II°°, H >4 b 

2 2 

we conclude that A = 1 and e = 1 + lg7r ^^+p • 


IV. CONFINEMENT 

The effective Lagrangian that we obtained by integrat¬ 
ing out the fermions has the form of a pure glue theory 
upto corrections ~ ^ and can be recast as 



with some simple rescalings of the fields and the coupling 
etc. The action of 18 will however have a speed of gluon 
given by 


1 _ 1 

y/Xe 

^ 3y/2nA 

gvi 


(19) 

( 20 ) 


as A = 1 and for fj,j A. Before elaborating on the 
scaling, we would like to point out that g', in eq. 18 
runs exactly like the coupling in a pure SU(3) Yang- 
Mills theory for momentums smaller than the gap. As 
the energy scales reach the gap from below and become 
bigger than the gap, the running behavior of the coupling 
will no longer be that of pure SU( 3) Yang-Mills and will 
get affected by the higher dimensional operators. Hence, 
the action of eq. 18 should be matched with the actions of 
eq. 13 at the scale A. To know how exactly this matching 
is done, we have to rescale time, fields and the coupling 
in eq. 13 as follows. We can rescale time direction, the 
gauge field and the coupling as t” } = At, A° a = v / eA g 
and g" = to rewrite the Lagrangian in eq. 13 as 


S 


1 

4 (g") 2 


d 4 x"(F") 2 


( 21 ) 


where a " = -C a s . Now we have to match the 

coupling constant in the action 18, a((g') with a" at a 
scale q' = A which in turn means 


«' S ( A ) = a 's 


a s {n) 


< a s . 


( 22 ) 


scale A to distinguish it from the zero density QCD con¬ 
finement scale Xqcd- We can figure out how the confine¬ 
ment scale of the effective theory compares with the gap 
by looking at the renormalization group equation 


A = A exp 


2 tt A 

bou'si A)J 


= AExp 


-2v^F- - - 

29Ay / o((A) 


(23) 


where b o at one loop with two flavours for color SU(3) 
is given by bo = 11 — ^rif = ^ and a' s is the coupling 
evaluated at the scale A. To evaluate the confinement 
scale we have to evaluate the gap. The dependence of A 
on the coupling is 

3tt 2 

A+ = b\fii\g(m) 5 exp(— 2g ^ ) ( 24 ) 

with b = 10 4 [1]. We compute and plot the confinement 
scale for a range of the isospin chemical potential from 
asymptotically large values where we trust our calcula¬ 
tions completely to regions of smaller values where the 
approximations start breaking down. In the region of 
smaller isospin density the confinement scale A could be 
appreciably different from what we found here. It would 
be instructive to compare our results with a future lattice 
calculation in the regime of small isospin density < 1000 
MeV where the different scales in the theory are not so 
well separated. Fig. 1 illustrates the various energy scales 
in the problem as a function of /i/. It also demonstrates 
how the scales get widely separated for asymptotically 
high isospin density where perturbation theory is valid 
and how they get closer in the opposite limit with low 
isospin density where the QCD coupling becomes con¬ 
siderably strong. For example if we look at fig. 1 for 
/ij ~ 5000 MeV, we notice that A ~ 400 MeV and this is 
an indication of the fact that our calculations are trust¬ 
worthy in this regime as the gap is much smaller than the 
isospin chemical potential. As expected from the eq. 23 
the confinement scale decreases exponentially with the 
isospin chemical potential and this is illustrated in fig. 2. 


V. EQUATION OF STATE AND THE CRITICAL 
TEMPERATURE OF DECONFINEMENT 
TRANSITION 


With this pure Yang-Mills theory in hand, we now con¬ 
centrate on the dynamics of the deconfinement transition. 
As mentioned in the introduction, a pure SU( 3) Yang- 
Mills theory undergoes a first order deconfinement phase 
transition with increasing temperature. This means that 
in our problem at finite temperature, for temperatures 
smaller compared to the gap there will be a deconfine¬ 
ment phase transition with increasing temperature. The 
energy scale of this deconfinement transition will be much 
smaller than the gap and we call this new confinement 


As mentioned in the introduction, it is possible to use 
our analysis to extract more information about the sys¬ 
tem than have already computed above. The key addi¬ 
tional information is the equation of state as a function 
of temperature in the regime of large /i/ and T <C A 
including the value of the critical temperature of the de- 
confinement transition. Of course, It is well known that 
even for pure Yang-Mills theory the equation of state and 
the critical temperature are not computable analytically 
due to the strength of the coupling. However, it can be 
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Hi in MeV 

FIG. 1: Confinement scale in MeV as a function of 
chemical potential in MeV 



Hi in MeV 


FIG. 2: The QCD confinement scale, the confinement 
scale at high isospin and the gap are shown in the same 
plot. The regime towards the right with higher isospin 
is where all these scales are well separated and we trust 
our calculations. The scales get closer as we go towards 
the left of the plot where the isospin chemical potential 
is small. This is the regime where a lattice calculation 
could determine how accurate the predictions are. 


computed on lattice using numerical techniques, and the 
lattice calculations to do are far simpler than for theories 
including dynamical quarks. If we parametrize the rela¬ 
tion between the transition for pure Yang-Mills with the 
confinement scale of the theory in the following way 

Tym = kymM-m (25) 

where Ty M is the critical temperature for the deconfine¬ 
ment transition, A ym is the confinement scale for pure 
Yang-Mills, kym , the proportionality constant between 
the two can be computed using lattice. The form of the 
relation between the critical temperature and confine¬ 
ment scale in Eq. 25 follows from dimensional analysis 
alone. 

We have already computed the confinement scale(A) 
for the low energy effective theory for our problem with 
QCD with two light quarks and a very high isospin chem¬ 
ical potential as a function of A qcd, the QCD confine¬ 
ment scale at zero density. Thus the critical temperature 


T c in our problem, will be related to A as 

T c = kym A. (26) 

The equation of state of QCD in the high isospin 
chemical potential limit with two light quark flavors and 
T <C A can also be related to that of pure Yang-Mills 
theory. In order to do this we need to discuss one ad¬ 
ditional aspect of the problem in detail. Other than the 
unscreened gluons at low energy scales there is also a 
Goldstone mode associated with the spontaneous break¬ 
ing of isospin rotations about the 3 direction in isospin 
space. The Goldstone mode corresponds to a rotation of 



The microscopic Lagrangian of Eq. 2 is invariant un¬ 
der this transformation but the condensate is not. 
Clearly this Goldstone mode contributes to the equation 
of state. However, it is also clear that the Goldstone 
mode does not interact with the dynamics of the gluons 
below the energy scale of the gap as the interactions are 
suppressed by powers of hi- Hence the deconfinement 
phase transition in the gluodynamics is not affected by 
the presence of this Goldstone mode. However, while 
the Goldstone mode does not interact with the glue- 
balls, it does contribute to the pressure as a function 
of temperature—that is the equation of state, and in a 
formal sense, the contribution from the Goldstone mode 
is of the same order in temperature as that of the gluons. 
Fortunately, as the gluons and the Goldstone mode do 
not interact and the two contributions can be computed 
separately and added to give the pressure. The contri¬ 
bution to the pressure from the Goldstone mode is given 

by 



where T is the temperature, v is the speed of the Gold¬ 
stone mode in the medium and is given by l/v3- This 
number l/\/3 is a result that is true for Goldstone modes 
in any finite density system [38]. The contribution com¬ 
ing from the gluons to the pressure is identical in form to 
that in pure Yang-Mills and can be computed on lattice. 
This contribution can be parametrized as 

fgluon = T 4 f (Vj (29) 

where / is a function of the temperature and the critical 
temperature. Here, note that the speed of gluons in the 
rescaled theory is given by 1. Again the form of the pres¬ 
sure in Eq. 29 can be obtained by dimensional analysis. 

The functional form of / (^ t T ''j as a function of T and 
Ty M can be computed using lattice. Hence we conclude 
that the equation of state for QCD with two light flavors 
of quark at a very high isospin chemical potential is given 
by 

/ T\ 7r 2 T 4 

p m = ryy + -^ ( 30) 







where / is a numerically calculable function of tempera¬ 
ture. 


VI. CONCLUSION 

The phenomenon of a deconfinement transition with 
increasing temperature at a very high isospin density 
might seem surprising in that a similar phenomenon is 
not expected at high baryon chemical potential (hb)- 
Given its unfamiliarity, it is useful to try to put in context 
what our results tell us about the nature of the phase dia¬ 
gram in the T — hi plane at zero /i b ■ Our calculations al¬ 
lowing the determination of the equation of state are only 
reliable in the high hi and low T region. Nevertheless, 
one can deduce a few qualitative features about the phase 
diagram. Several features of the phase diagram should 
be regarded as being on solid ground. The first is that at 
Hi = 0, there are no phase transitions as T changes. The 
second is that the for sufficiently high hi there exists a 
superconducting phase with a gap serving as an order pa¬ 
rameter. At low temperatures and isospin densities this 
superconducting phase may be regarded as a pion con¬ 
densate. At high temperature however this condensate 
melts. Thus there must exist a curve which separates 
the superconducting phase from the normal phase. The 
third robust feature of the phase diagram was the fo¬ 
cus of this paper—namely the existence of a first-order 
phase transition at sufficiently high hi between a con¬ 
fining and a deconfining phase. This implies the exis¬ 
tence of a curve which separates the confined phase from 
the deconfined phase. The critical temperature of the 
deconfinement transition is proportional to the confine¬ 
ment scale we have evaluated above as written in Eq. 26. 
As was seen in fig. 3c it decreases exponentially with 
increasing isospin density in the regime where we trust 
our calculations. But in the regime where our calcula¬ 
tions are not valid the behavior of this line is a matter of 
speculation. 

These three facts imply the existence of certain spe¬ 
cial points in the phase diagram. For the moment let 
us assume that regardless of /i/, there are at most two 
phase transitions as a function of T, one associated with 
superconductivity and the other associated with confine¬ 
ment. However, as was pointed out in [1] we also know 
that curves associated with these transitions must end 
somewhere at sufficiently low /i/, since there are no tran¬ 
sitions at hi = 0. The points where the curves end are 
special. For the case of the superconducting transition, 
this special point occurs at T = 0 and hi = m* as is 
seen in all three curves in Fig. 3. An interesting and 
nontrivial fact is that one can infer the existence of an 
additional special point where the deconfinement tran¬ 
sition ends. Son and Stephanov speculated about the 
nature of this point [1]. They argued that it is very plau¬ 
sible that there is no phase transition at zero T as one 
increases hi from m n (where the superconducting transi¬ 
tion occurs) as there is no symmetry. A very natural way 


T 





FIG. 3: The three possibile scenarios describing what 
happens to the deconfinement phase transition line as 
Hi is lowered. T^ c corresponds to the critical 
temperature for the superconducting transition and T c 
is the critical temperature for the deconhnement 
transition. 


for this to occur is the one suggested in [1] and shown in 
Fig. 3a, namely for the curve of first-order transitions 
separating the confined from deconfined phases to end 
at a critical point. This possibility is particularly ex¬ 
citing since critical points are characterized by universal 
behavior. However, we should note that it is not guaran- 
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teed to occur. For example, one could imagine that the 
special point associated with the end of the confinement 
transition could be a triple point as shown in Fig. 3b 
Alternatively, one might imagine, that despite the plau¬ 
sibility argument in ref. [1], that the curve might end 
at T = 0 as in Fig. 3c indicating a second transition 
at T = 0. One can also imagine more baroque scenar¬ 
ios with additional phases. The point, though is that 
however, the curve ends, it must end somewhere and the 
point at which this happens is intrinsically of interest. 
Ultimately, one would hope that lattice calculations can 


determine the phase structure and the position of these 
interesting points. 
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